Topological Charge Fluctuations and Low-Lying Dirac Eigenmodes 
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We discuss the utility of low-lying Dirac eigenmodes for studying the nature of topological charge fluctuations 
in QCD. The implications of previous results using the local chirality histogram method are discussed, and 
the new results using the overlap Dirac operator in Wilson gauge backgrounds at lattice spacings ranging from 
a « 0.04 fm to a f» 0.12 fm are reported. While the degree of local chirality does not change appreciably closer to 
the continuum limit, we find that the size and density of local structures responsible for chiral peaking do change 
significantly. The resulting values are in disagreement with the assumptions of the Instanton Liquid Model. We 
conclude that the fluctuations of topological charge in the QCD vacuum are not locally quantized. 



The instanton (semiclassical) picture of the 
QCD vacuum has undergone several refinements 
since its beginnings as a dilute instanton gas, 
but there are two cornerstones that are inherent 
in this approach. First, it is imagined that the 
gauge field fluctuates in such a way as to form 
the lumps of approximately quantized topologi- 
cal charge. More precisely, it is assumed that 
there are regions of strong gauge fields that can 
be enclosed by a hypersurface on which the gauge 
field is approximately pure gauge. The topologi- 
cal charge contained in such regions thus have ap- 
proximately quantized values, and we will refer to 
this situation as "lumpy vacuum" or vacuum with 
locally quantized topological charge. Secondly, 
the gauge field in such lumps is assumed to be 
overwhelmingly (anti)self-dual. Rather success- 
ful phenomenology, the Instanton Liquid Model 
(ILM) 1^, has been built around these assump- 
tions. The basic characteristics of the ILM are the 
average instanton radius p « 1/3 fm and density 
n « 1 fm"'' with values rather tightly restricted. 

Using fermionic response as a signal for gauge 
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fluctuations is a plausible strategy on the lattice, 
since lattice gauge fields are quite rough and be- 
lieved to be affected strongly by artifacts at the 
scale of the lattice spacing. The attemps to elim- 
inate the unphysical part of UV fluctuations di- 
rectly from the gauge fleld are quite subjective 
and significantly biased as to the local structure 
of the vacuum they reveal. Indeed, cooling in- 
evitably leads to lumpy and self-dual gauge fields 
if global topology is preserved in the process. 
Similarly, smoothing is biased towards the lumpy 
vacuum but not necessarily towards self-duality. 

The eigenset of the Dirac operator fully en- 
codes the information about the topological fluc- 
tuations of the underlying gauge field. While re- 
lying on the complete eigenset is impractical, the 
importance of these fluctuations mainly rests in 
their influence on the propagation of light quarks 
which, in turn, depends crucially on the low-lying 
modes. Hence, it is plausible to look for the im- 
prints of topological charge fluctuations in Dirac 
near-zero modes This approach actually rep- 
resents a form of smoothing because the infrared 
modes are least sensitive to the short-distance fea- 
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tures atop larger structures in the gauge poten- 
tial. However, "fermionic smoothing" is deter- 
mined by the relevant fermionic dynamics, thus 
eliminating the subjective element. 

While lattice fermion is not unique (different 
fermions smooth differently), all good actions 
should give consistent results sufhciently close to 
the continuum limit. However, to make meaning- 
ful comparisons, the same fermionic action should 
be used as the lattice spacing changes. There ap- 
pears to be an interesting complementarity here 
of the type "degree of chiral symmetry x reso- 
lution = consf . Valid ultralocal actions can not 
have exact chiral symmetry ||] , but they can have 
perfect lattice "resolution" in the sense that the 
fermion only feels the gauge potential at its lo- 
cation. The actions with exact chiral symmetry, 
however, must be non-ultralocal and their reso- 
lution is poor. The advantage of the overlap op- 
erator over the Wilson-Dirac operator which still 
reflects the global topology quite well and has per- 
fect resolution, is thus not obvious a ■priori. 

A concrete proposal along these lines is the in- 
vestigation of chiral orientation parameter X{n) 
in the low eigenmodes X{n) is a local an- 
gle in the IV'lI — \i'R\ plane, normalized to lie in 
the interval [—1, 1]. This is quite useful because 
one can argue that the distribution of this quan- 
tity over active points on the lattice will have a 
predictable qualitative behaviour when the field 
is lumpy or self-dual. For a lumpy configuration, 
the emergence of topological near- zero-modes is 
expected, and the distribution of X should be 
strongly peaked at the extremal values. When 
the gauge field exhibits extended regions of strong 
self-dual fields, the distribution should tend to 
be peaked as well. We emphasize that these two 
cases can happen independently^ e.g. the lumpy 
field does not have to be self-dual for the dis- 
tribution to be peaked. Consequently, such a 
qualitative behaviour is a necessary condition for 
instanton-like vacuum but by no means sufficient! 

The first study of local chirality with Wilson 
fermions on Wilson gauge backgrounds at lat- 
tice spacing aw 0.17 fm found results inconsistent 
with the instanton-like vacuum ||^. Later works 
with Wilson fermions at smaller lattice spacing 
and also with the overlap revealed however, that 



a visible degree of double-peaked behaviour is 
present Our overlap calculation on Wilson 
gauge fields at a ~ 0.123, 0.042 fm gave the fol- 
lowing results 



Chiral Orientation Parameter Distribution 
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indicating that there does not appear to be a 
significant change in the degree of local chirality 
closer to the continuum limit. While the observed 
peaking does not imply the dominance of instan- 
tons, it does force us to look for quantitative com- 
parisons rather than qualitative ones. It might be 
indeed quite revealing to compare the probability 
distribution of X for an instanton liquid ensemble 
and for QCD using the overlap operator. To il- 
lustrate this point, we compare our results to the 
prediction of the dilute instanton gas model, 

beta=6.55 versus Instanton Gas 



I 0-2 




confirming that it is not applicable. 

If the low-lying mode arises from mixing of 
t'Hooft "would be" zero modes associated with 
instantons, then the regions of coherent local chi- 
rality should be of a particular size and abun- 
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dance. To check whether the behaviour of X- 
distribution can be ascribed to gauge structures 
with ILM parameters, we have found the focal 
maxima over the distance V3 of density {ip~^ilj)i 
and sefocted ones that (1) decay on average afong 
the eight axis directions through distance -^3 
(this excludes accidental maxima and those not 
likely due to instantons), and (2) contribute to 
the peaks of X-distribution. We then determined 
the size of selected structures as (a) the radius of 
the largest hypersphere centered at the maximum 
and containing the sites with the same sign of chi- 
rality {i/j'^ j^tp) , and (b) the radius of the largest 
hypersphere within which density doesn't fall be- 
low 1/8 of the center value (this corresponds to 
the radius in the case of an instanton). Using the 
first nonzero modes we found 

Size from Chirality Sign and Density Decay 




with values significantly below the ILM assump- 
tion close to continuum limit. Similarly, for the 
density of these structures we have obtained 
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in a marked disagreement with the ILM. A 
few points need to be emphasized here. (A) 
There should not be any additional smoothing or 
other manipulation done on the low-lying eigen- 
modes. The lowest nonzero infrared mode pro- 
vides a maximally smoothed imprint of topolog- 
ical charge fluctuations as seen by the propagat- 
ing fcrmion. (B) The typical "size" of the chiral- 
ity correlator in the lowest modes (representing 
an upper bound on the average size of underly- 
ing individual structures) at aw 0.042 fm is about 
0.23 fm. (C) The possibility that the large den- 
sity of structures we find is due to dislocations is 
unlikely and can be excluded in several ways. 

How do we interpret our data? They are 
easy to understand if one abandons the assump- 
tion that topological charge in QCD vacuum is 
approximately locally quantized. For quantized 
topological charges, the calculable occurence of 

0. 15 fm instantons should be negligible, but there 
is no such restriction for inhomogeneous fluctu- 
ations without quantization. Similarly, if the 
large number of structures we found corresponded 
to quantized "lumps" , then the typical values 
of global topological charge (and susceptibility) 
would be much higher than we observe, but not 
necessarily so without quantization. We conclude 
that it is unlikely that the observed peaking of lo- 
cal chirality distribution is due to instanton-like 
objects with the parameters of the ILM. Instead, 
we find an evidence for serious breakdown of the 
semiclassical picture of the QCD vacuum Q. 
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